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Abstract 

The Schrodinger cat male and female states are discussed. The Wigner and Q-functions 
of generalized correlated light are given. Linear transformator of photon statistics is reviewed. 


1 Introduction 

The integral of motion which is quadratic in position and momentum was found for classical 
oscillator with time-dependent frequency by Ermakov [1]. Two time-dependent integrals of motion 
which axe linear forms in position and momentum for the classical and quantum oscillator with 
time-dependent frequency were found in [2]; for a charge moving in varying in time uniform 
magnetic field, this was done in [3]. For the multimode nonstationary oscillatory systems, such 
new integrals of motion, both of Ermakov’s type (quadratic in positions and momenta) and linear 
in position and momenta, generalizing the results of [2] were constructed in [4]. We will consider 
below the parametric oscillator using the integrals of motion.The Wigner function of multimode 
squeezed light is studied using such special functions as multivariable Hermite polynomials. The 
theory of parametric os cillator is appropriate to consider the problem of creation of photons 
from vacuum in a resonator with moving walls (with moving mirrors) which is the phenomenon 
based on the existence of Casimir forces (so-called nonstationary Casimir effect). The resonator 
with moving boundaries (moving mirrors, media with time-dependent refractive index) produces 
also effect of squeezing in the light quadratures. In the high energy physics very fast particle 
collisions may produce new types of states of boson fields (pions, for example) which are squeezed 
and correlated states studied in quantum optics but almost unknown in particle physics, both 
theoretically and experimentally. 


2 Multimode Quadratic Systems 

The generic nonstationary linear system has the Hamiltonian 

H = ^QB{t)Q + C(t)Q, (1) 

where we use 2N-vectors Q = (pi, P2> ■ • • iPn> 9 i > Q2> • ■ ■ >Qn) an d C(t), as well as 2Nx2N- 
matrix B{t), the Planck constant K = 1. This system has 2N linear integrals of motion [5], [6] 
which may be written in vector form 

Qo(«) = A(t)Q + A (t). (2) 
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The real symplectic matrix A (t) is the solution to the system of equations 

A (t) = A(t)EB(t), 

A(0) = 1, (3) 

where the real antisymmetrical matrix E is 2N-dimensional analog of the Pauli matrix ia y , and 
the vector A (t) is the solution to the system of equations 

A (t) - A(0EC(f), 

A(0) = 0. (4) 

If for time t = 0, one has the initial Wigner function of the system in the form 

W{p,q,t = 0) = Wo(Q), ( 5 ) 

the Wigner function of the system at time t is (due to the density operator is the integral of 
motion) 

W(p,q,«) = W 0 [A(*)Q + A(*)]. (6) 

This formula may be interpreted as transformation of input Wigner function into output Wigner 
function due to symplectic quadrature transform (2). An optical linear transformator of photon 
distribution function using this output Wigner function i s suggested in [7]. 

The Hamiltonian (1) may be rewritten in terms of creation and annihilation operators 

H = \AD(t)A + E(t)A, (7) 

where we use 2N-vectors A = ( Oi, 02 , . . . , aw, ait, ° 2 t, • • ■ , <Pvt) and E (t), as well as 2Nx2N- 
matrix D(t). This system has 2N linear integrals of motion [5], [6] which are written in vector 
form 

Ac (t) = M{t) A + N (t). (8) 

The complex matrix M(t) is the solution to the system of equations 

M(t) = M(t)aD(t), 

M(0) = 1, (9) 

where the imaginary antisymmetric matrix o is 2Nx2N-analog of the Pauli matrix — <r v , and 

the vector N(t) is the solution to the system of equations 

N(f) = M{t)aE(t), 

N(0) = 0. (10) 

Analogously to the Wigner function evolution, if for time t = 0, one has the initial Q-function 
of the system in the form 


Q(a,a*,t = 0) = Q 0 (A), A=(Q! ) Q! , '), 


( 11 ) 
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( 12 ) 


the Q-function of the system at time t is 

Q(a,a\t) = Q 0 [M(t)A + m)]- 

Here a = ( q + ip)/\/2. 

For time-independent Hamiltonian (1), the matrix A (t) is 

A(f) = exp(EBt), (13) 

and the vector A (t) is 

A (t) = /‘exp(E Bt) E C(r) dr. (14) 

Jo 

For time-independent Hamiltonian (7), the matrix M(t ) is 

M(t) = exp(cr-Dt), (15) 


and the vector N(£) is 

N(f) = / exp (a Dr) a E(r) dr. (16) 

Jo 

For time-dependent linear systems, the Wigner function of generic squeezed and correlated state 
(generalized correlated state [8] ) has Gaussian form and it was calculated in [5]. 

Thus the evolution of the Wigner function and Q-function for systems with quadratic Hamilto- 
nians for any state is given by the following prescription. Given the Wigner function W(p, q, t = 0) 
for the initial moment of time t = 0. Then the Wigner function for time t is obtained by the 


replacement 


W(p,q,t)= W(p(t), q(t), t = 0), 


where time-dependent arguments are the linear integrals of motion of the quadratic system found 
in [5], [4], and [9]. This formula was given as integral with 6-function kernel in [10]. The linear 
integrals of motion describe initial values of classical trajectories in the phase space of the system. 
The same ansatz is used for the Q-function. Namely, given the Q-function of the quadratic system 
Q(B, t = 0) for the initial moment of time t = 0. Then the Q-function for time t is given by 


the replacement 


Q(B, t)= i = 0), 


where the 2N-vector B (t) is the integral of motion linear in the annihilation and creation op- 
erators. This ansatz follows from the statement that the density operator of the Hamiltonian 
system is the integral of motion, and its matrix elements in any basis must depend on appropriate 
integrals of motion. 


3 Multimode Mixed Correlated Light 

The most general mixed squeezed state of the N-mode light with a Gaussian density operator p 
is described by the Wigner function VF(p,q) of the generic Gaussian form, 

W(p,q) - detMexp j-^(Q- < Q >)M _1 (Q~ < Q >) > ( 17 ) 
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where 2N parameters < Pi > and < qi >, i = 1,2, . . . ,N, combined into vector < Q>, are 
average values of quadratures, 

< p > = TV pp, 

< q > = TV pq. (18) 

A real symmetric dispersion matrix M consists of 2N 2 +N variances 

Ma0 = \ < QaQp + QpQa > - < Qa >< Qff >, a, 0 = 1 , 2, . . . , 2N. (19) 

They obey uncertainty relations constraints [5]. According to previous section the Wigner function 
of parametric linear system with initial value (17) is 

W (p> t) = det M exp -^(A(t)Q + A(t)~ < Q >)M -1 (A(*)Q 4- A(t)- < Q >) , (20) 

The photon distribution function of the state (17) 


V n = TV p|n >< n|, n = (m,n 2 , . . . .n^), (21) 

where the state |n > is photon number state, which was calculated in [11], [12] and it is 

#nn } (v) 

Pn-Po— (22) 

The trace (21) may be calculated using the explicit form of the Wigner function of the opera- 
tor |m >< n| (see, [5]) which is the product of Wigner functions of one-dimensional oscillator 
expressed in terms of Laguerre polynomials of the form 


W mn (p,q) = 2 m - n+l (-l) n 


n\ ( q — ip\ m ” 


ml V v/2 


e -tfW) L rn-n ( 2 (g 2 + p 2 )) . 


rni 

The function H A n J (y) is multidimensional Hermite polynomial. The probability to have no 
photons is 

T ( 1 M _1/2 r ~\ 

Vo = det + -I 2 nJ exp [- < Q > (2M + 1 2 n) 1 < Q > , (24) 

where we introduced the matrix 

R- = 2U^(1 + 2M) _1 U* — (Tn x , (25) 

and the matrix 

° N * = ( I* *0 ) ' (26) 

The argument of Hermite polynomial is 


y = 2U*(I !W - 2M)" 1 < Q >, 


(27) 
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and the 2N-dimensional unitary matrix 


v= \(-a N a„\ ( 28 ) 

■y/2 \ lyv I N ) 

is introduced, in which Iyv is the NxN— identity matrix. Also, we use the notation 

n! = Tiling! • • • n^!. 

The mean photon number for j-th mode is expressed in terms of photon quadrature means 
and dispersions 

< Tlj >= -{<7 PjPj + <?q jqj ~ ^) + ^(< Pj > 2 + < Qj >2 )- (29) 

The photon distribution function for transformed state (20) is given by the same formulae (22), 
(24)-(28) but with changed dispersion matrix 

M = A -1 MA _lt , (30) 


and quadrature means 

< Q >= A -1 (A— < Q >). (31) 

Thus we have a linear transformator of photon statistics suggested in [7]. 

Let us now introduce a complex 2N-vector B = (/?i, pi, ■ Pn, Pit Pi, •••> Pn)- Then 
the Q-function is the diagonal matrix element of the density operator in coherent state basis 
\ Pit Pit ■■■ t Pn > ■ This function is the generating function for matrix elements of the density 
operator in the Fock basis |n> which has been calculated in [12]. In notations corresponding to 
the Wigner function (17) the Q-function is 


Q(B) = Vo exp 


-]-B(R + a Nx )B + BRy . 
z 


(32) 


Thus, if the Wigner function (17) is given one has the Q-function. Also, if one has the Q-function 
(32), i.e., the matrix R and vector y, the Wigner function may be obtained due to relations 

M = U*(fl + ajv*) _1 ut-l/2, 

< Q > = U*[l - (fl + <7jvx)" V N x]y. (33) 


Multivariable Hermite polynomials describe the photon distribution function for the multimode 
mixed and pure correlated light [11], [13], [14]. The nonclassical state of light may be created due 
to nonstationary Casimir effect [15], [16] and the multimode oscillator is the model to describe 
the behaviour of squeezed and correlated photons. 


4 Parametric Oscillator 


For the parametric oscillator with the Hamiltonian 


d 2 uPjt)x^ 

2dx 2 2 ’ 


(34) 
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where we take h = m = w(0) = 1, there exists the time-dependent integral of motion found in [2] 


A = 


V2 


[e(t)p - e{t)x]. 


(35) 


where 

e(t) + u> 2 (t)e(t) = 0, e(0) = 1, e(0) = z, (36) 

satisfying the commutation relation 

[A, A\\ = 1. (37) 

It is easy to show that packet solutions of the Schrodinger equation may be introduced and 
interpreted as coherent states [2], since they are eigenstates of the operator A (35), of the form 


ty Q (x,t) = tyo(x,t)exp 



a 2 e*(t ) y/2ax 
2 e(t) + e 


(38) 


where 


'I'oOM) 


= 7r- 1/4 e(t)- 1/2 exp 


ie(t)x 2 

2e(0 


(39) 


is analog of the ground state of the oscillator and a is a complex number. 

Variances of the position and momentum of the parametric oscillator in the state (38), (39) 


are 



(40) 


and the correlation coefficient r of the position and momentum has the value corresponding to 
minimization of the Schrodinger uncertainty relation [17] 


O xOp — 


l l 
41 -r 2 ' 


(41) 


If o x < 1/2 (<r p < 1/2) we have squeezing in photon quadrature components. 

The analogs of orthogonal and complete system of states which are excited states of stationary 
oscillator are obtained by expansion of (38) into power series in a. We have 

= ( 2 ^) (roi) ’ (42) 

and these squeezed and correlated number states are eigenstates of invariant A. In case of 
periodical dependence of frequency on time the classical solution in stable regime may be taken 
in Floquet form 

e(t) = e^u(t), (43) 

where u(t) is a periodical function of time. Then the states (42) are quasienergy states realizing 
the unitary irreducible representation of time translation symmetry group of the Hamiltonian 
and the parameter k determines the quasicnergy spectrum. Unstable classical solutions give 
continuous spectrum of quasienergy states. 
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The partial cases of parametric oscillator are free motion ( ui(t) — 0 ), stationary harmonic 
oscillator ( u 2 (t) = 1 ), and repulsive oscillator ( cn 2 (t) = -1 ). The solutions obtained above are 
described by the function e(t) which is equal to e(t) = 1 + it, for free particle, e(t) = e l£ , for 
usual oscillator, and e(t) = cosht + rsinhf, for repulsive oscillator. 

Another normalized solution to the Schrodinger equation 


am (*£ > 1} 


2N m 'b 0 {x,t) exp 



e^tW\ 
2e{t) ) 


cosh 


y/2 ax 

~W’ 


(44) 


where 


exp(|a| 2 /2) 

2ycosh |a| 2 


(45) 


is the even coherent state [18] (the Schrodinger cat male state). The odd coherent state of the 
parametric oscillator (the Schrodinger cat female state) 


S' Q f(x,t ) = 2N f ty 0 (x,t)exp 



e*(t)a 2 \ 

2e(t) ) 


sinh 


\[2a.x 


(46) 


where 


N f = 


exp(|a| 2 /2) 
2^/sinh |a| 2 


(47) 


satisfies the Schrodinger equation and is the eigenstate of the integral of motion A 2 (as well as the 
even coherent state) with the eigenvalue a 2 . These states are one-mode examples of squeezed and 
correlated Schrodinger cat states constructed in [19]. The experimental creation of the Schrodinger 
cat states is discussed in [20]. These states belo ng to family of nonclassical superposition states 
studied in [21], [22]. 
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